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Abstract 

We address the connectivity of large-scale ad hoc heterogeneous wireless networks, where secondary 
users exploit channels temporarily unused by primary users and the existence of a communication link 
between two secondary users depends on not only the distance between them but also the transmitting 
and receiving activities of nearby primary users. We introduce the concept of connectivity region defined 
as the set of density pairs — the density of secondary users and the density of primary transmitters 
— under which the secondary network is connected. Using theories and techniques from continuum 
percolation, we analytically characterize the connectivity region of the secondary network and reveal 
the tradeoff between proximity (the number of neighbors) and the occurrence of spectrum opportunities. 
Specifically, we establish three basic properties of the connectivity region - contiguity, monotonicity of 
the boundary, and uniqueness of the infinite connected component, where the uniqueness implies the 
occurrence of a phase transition phenomenon in terms of the almost sure existence of either zero or 
one infinite connected component; we identify and analyze two critical densities which jointly specify 
the profile as well as an outer bound on the connectivity region; we study the impacts of secondary 
users' transmission power on the connectivity region and the conditional average degree of a secondary 
user, and demonstrate that matching the interference ranges of the primary and the secondary networks 
maximizes the tolerance of the secondary network to the primary traffic load. Furthermore, we establish 
a necessary condition and a sufficient condition for connectivity, which lead to an outer bound and an 
inner bound on the connectivity region. 
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I. Introduction 

The communication infrastructure is becoming increasingly heterogeneous, with a dynamic 
composition of interdependent, interactive, and hierarchical network components with different 
priorities and service requirements. One example is the cognitive radio technology dll] for 
opportunistic spectrum access which adopts a hierarchical structure for resource sharing [Q]- 
Specifically, a secondary network is overlaid with a primary network, where secondary users 



identify and exploit temporari 



y and locally unused channels without causing unacceptable 



interference to primary users [|2|]. 

A. Connectivity and Connectivity Region 

While the connectivity of homogeneous ad hoc networks consisting of peer users has been 

nTinnnnuri 

well studied (see, for example, DSL |4|, |5|, la 13, ISL |9|, 110|]), little is known about the connectivity 



of heterogeneous networks. The problem is fundamentally different from its counterpart in 
homogeneous networks. In particular, the connectivity of the low -priority network component 
depends on the characteristics (traffic pattern/load, topology, interference tolerance, etc.) of the 
high-priority component, thus creating a much more diverse and complex design space. 

Using theories and techniques from continuum percolation, we analytically characterize the 
connectivity of the secondary network in a large-scale ad hoc heterogeneous network. Specifi- 
cally, we consider a Poisson distributed secondary network overlaid with a Poisson distributed 
primary network in an infinite two-dimensional Euclidean spacqlJ. We define network connectivity 
as the existence of an infinite connected component almost surely (a.s.), i.e., the occurrence of 
percolation. We say that the secondary network is strongly connected when it contains a unique 
infinite connected component a.s. 

Due to the hierarchical structure of spectrum sharing, a communication link exists between 
two secondary users if the following two conditions hold: (CI) they are within each other's 
transmission range; (C2) they see a spectrum opportunity determined by the transmitting and 
receiving activities of nearby primary users (see Sec. III-Bll) . Thus, given the transmission power 

'This infinite network model is equivalent in distribution to the limit of a sequence of finite networks with a fixed density 
as the area of the network increases to infinity, i.e., the so-called extended network \\\ . It follows from the arguments similar 
to the ones used in |l^ Chapter 3] for homogeneous ad hoc networks that this infinite ad hoc heterogeneous network model 
represents the limiting behavior of large-scale networks. 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY, AUGUST, 2009. 



3 



and the interference tolerance of both the primary and the secondary users, the connectivity of 
the secondary network depends on the density of secondary users (due to (CI)) and the traffic 
load of primary users (due to (C2)). 

We thus introduce the concept of connectivity region C, defined as the set of density pairs 
(As, Xpt) under which the secondary network is connected, where Ag denotes the density of the 
secondary users and Apr the density of primary transmitters (representing the traffic load of the 
primary users). As illustrated in Fig. [H a secondary network with a density pair (A5, Apr) inside 
this region is connected: the secondary network has a giant connected component which includes 
infinite secondary users. The existence of the giant connected component enables bidirectional 
communications between distant secondary users via multihop relaying. On the other hand, a 
secondary network with a density pair (A^, Apt) outside this region is not connected: the network 
is separated into an infinite number of finite connected components. Consequently, any secondary 
user can only communicate with users within a limited range. 



PT 




A^ \s 

Fig. 1. The connectivity region C (the upper boundary Apy(As) is defined as the supremum density of the primary transmitters 
to ensure connectivity with & fixed density of the secondary users; the critical density of the secondary users is defined as the 
infimum density of the secondary users to ensure connectivity under a positive density of the primary transmitters; the critical 
density Apj, of the primary transmitters the supremum density of the primary transmitters to ensure connectivity with & finite 
density of the secondary users). 

The objective of this paper is to establish analytical characterizations of the connectivity region 
and to study the impact of system design parameters (in particular, the transmission power of the 
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secondary users) on the network connectivity. Main results are summarized in the subsequent 
two subsections. 

B. Analytical Characterizations of the Connectivity Region 

We first establish three basic properties of the connectivity region: contiguity, monotonicity 
of the boundary, and uniqueness of the infinite connected component. Specifically, based on a 
coupling argument, we show that the connectivity region is a contiguous area bounded below 
by the A5-axis and bounded above by a monotonically increasing function \*pt{^s) (see Fig. [T]), 
where the upper boundary Apj^(As') is defined as 

\*prp[\s) = supjApT : G{^s, ^pt) is connected.}, 

with G{\s, Apt) denoting the secondary network of density A5 overlaid with a primary network 
specified by the density Ap^ of the primary transmitters. The uniqueness of the infinite connected 
component is established based on the ergodic theory and certain combinatorial results. It shows 
that once the secondary network is connected, it is strongly connected. 

Second, we identify and analyze two critical parameters of the connectivity region: A^ and 



Ap^. They jointly specify the profile as well as an outer bound on the connectivity region. 
Referred to as the critical density of the secondary users, A^ is the infimum density of the 
secondary users to ensure connectivity under a positive density of the primary transmitters: 

A5= infjA^ : 3ApT > s.t. ^(A5,Apr) is connected}. 

We show that A J equals the critical density Ac of a homogeneous ad hoc network (i.e., in the 



absence of primary users), which has been well studied 111311 . This result shows that the "takeoff" 
point in the connectivity region is completely determined by the effect of proximity — the number 
of neighbors (nodes within the transmission range of a secondary user). 



Referred to as the critical density of the primary transmitters, Ap-^ is the supremum density 
of the primary transmitters to ensure the connectivity of the secondary network with a finite 
density of the secondary users: 



X*PT= sup{Apr : 3As < 00 s.t. ^(As,Apt) is connected}. 



We obtain an upper bound on Ap^ which is shown to be achievable in simulations. More 
importantly, this result shows that when the density of the primary transmitters is higher than 
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the (finite) value given by this upper bound, the secondary network cannot be connected no matter 
how dense it is. This parameter Ap^ thus characterizes the impact of opportunity occurrence on 
the connectivity of the secondary network: when the density of the primary transmitters is beyond 
a certain level, there are simply not enough spectrum opportunities for any secondary network 
to be connected. 

Since a precise characterization of the upper boundary \*prp{Xs) of the connectivity region is 
intractable, we establish a necessary and a sufficient condition for connectivity to provide an 
outer and an inner bound on the connectivity region. The necessary condition is expressed in the 
form of the conditional average degree of a secondary user, and is derived by the construction of 
a branching process. The sufficient condition is obtained by the discretization of the continuum 
percolation model into a dependent site percolation model. 




Fig. 2. Simulated connectivity regions for two different transmission powers {ptx denotes the transmission power of the 
secondary users, and the large ptx is 3" times the small ptx, where a is the path-loss exponent). 

C. Impact of Transmission Power on Connectivity: Proximity vs. Opportunity 

The study on the impact of the secondary users' transmission power on the network con- 
nectivity reveals an interesting tradeoff between proximity and opportunity in the design of 
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heterogeneous networks. As illustrated in Fig. [2l we show that increasing ptj. enlarges the 
connectivity region C in the A^-axis {i.e., better proximity leads to a smaller "takeoff" point), 
but at the price of reducing C in the Ap^-axis. Specifically, with a large p^^., few secondary 
users experience spectrum opportunities due to their large interference range with respect to the 
primary users. This leads to a poor tolerance to the primary traffic load parameterized by XpT- 
The transmission power ptx of the secondary network should thus be chosen according to 
the operating point of the heterogeneous network given by the density of the secondary users 
and the traffic load of the co-existing primary users. Using the tolerance to the primary traffic 
load as the performance measure, we show that the interference range r/ of the secondary users 
should be equal to the interference range Ri of the primary users in order to maximize the upper 
bound on the critical density \*prp of the primary transmitters. Given the interference tolerance 
of the primary and secondary users, we can then design the optimal transmission power ptx of 
the secondary users based on that of the primary users. 



D. Related Work 

To our best knowledge, the connectivity of large-scale ad hoc heterogeneous networks has 
not been characterized analytically or experimentally in the literature. There are a number of 
classic results on the connectivity of homogeneous ad hoc networks. For example, it has been 
shown that to ensure either 1 -connectivity (there exists a path between any pair of nodes)J5, ^ 
or /^-connectivity (there exist at least k node-disjoint paths between any pair of nodes) ||8|], the 
average number of neighbors of each node must increase with the network size. On the other 
hand, to maintain a weaker connectivity - p-connectivity {i.e., the probability that any pair of 
nodes is connected is at least p), the average number of neighbors is only required to be above 
a certain 'magic number' which does not depend on the network size 

The theory of continuum percolation has been used by Dousse et al. in analyzing the con- 
nectivity of a homogeneous ad hoc network under the worst case mutual interference yl, |4i]. 



In m, 



1011 . the connectivity and the transmission delay in a homogeneous ad hoc network with 



statically or dynamically on-off links are investigated from a percolation-based perspective. 

The optimal power control in heterogeneous networks has been studied in [|l4ll . which focuses 
on a single pair of secondary users in a Poisson network of primary users. The impacts of sec- 
ondary users' transmission power on the occurrence of spectrum opportunities and the reliability 
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of opportunity detection are analytically characterized. 

E. Organization and Notations 

The rest of this paper is organized as follows. Sec. |ll] presents the Poisson model of the 
heterogeneous network. In particular, the conditions for the existence a communication link in 
the secondary network is specified based on a rigorous definition of spectrum opportunity. In 
Sec. Uni we introduce the concept of connectivity region and establish its three basic properties. 
The two critical densities are analyzed, followed by a necessary and a sufficient condition for 
connectivity. In Sec. |IVl we demonstrate the tradeoff between proximity and opportunity by 
studying the impacts of the secondary users' transmission power on the connectivity region and 
on the conditional degree of a secondary user. The optimal transmission power of the secondary 
users is obtained under the performance measure of the secondary network's tolerance to the 
primary traffic load. Sec |V] contains the detailed proofs of the main results, and Sec. IVTl concludes 
the paper. 

Throughout the paper, we use capital letters for parameters of the primary users and lowercase 
letters for the secondary users. 

II. Netvv'Ork Model 

We consider a Poisson distributed secondary network overlaid with a Poisson distributed 
primary network in an infinite two-dimensional Euclidean space. The models of the primary and 
secondary networks are specified in the following two subsections. 

A. The Primary Network 

The primary transmitters are distributed according to a two-dimensional Poisson point process 
with density Apt. To each primary transmitter, its receiver is uniformly distributed within 
its transmission range Rp. Here we have assumed that all primary transmitters use the same 
transmission power and the transmitted signals undergo an isotropic path loss. Based on the 
displacement theorem [|l5l . Chapter 5], it is easy to see that the primary receivers form a two- 
dimensional Poisson point process with density Apt. Note that the two Poisson processes formed 
by the primary transmitters and receivers are correlated. 
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B. The Secondary Network 

The secondary users are distributed according to a two-dimensional Poisson point process 
with density A5, independent of the Poisson processes of the primary transmitters and receivers. 
The transmission range of the secondary users is denoted by Vp. 

1 ) Communication Links: In contrast to the case in a homogeneous network, the existence of 
a communication link between two secondary users depends on not only the distance between 
them but also the availability of the communication channel {i.e., the presence of a spectrum 
opportunity). The latter is determined by the transmitting and receiving activities in the primary 
network as described below. 

As illustrated in Fig.[3l there exists an opportunity from A, the secondary transmitter, to B, the 
secondary receiver, if the transmission from A does not interfere with nearby primary receivers 
in the solid circle, and the reception at B is not affected by nearby primary transmitters in the 
dashed circle jl^. Referred to as the interference range of the secondary users, the radius vj 
of the solid circle at A depends on the transmission power of A and the interference tolerance 
of the primary receivers, whereas the radius Rj of the dashed circle (the interference range of 
the primary users) depends on the transmission power of the primary users and the interference 
tolerance of B. 




Fig. 3. Definition of spectrum opportunity. 



It is clear from the above discussion that spectrum opportunities depend on both transmitting 
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and receiving activities of the primary users. Furthermore, spectrum opportunities are asymmetric. 
Specifically, a channel that is an opportunity when A is the transmitter and B the receiver may 
not be an opportunity when B is the transmitter and A the receiver. In other words, there 
exist unidirectional communication links in the secondary network. Since unidirectional links 



are difficult to utilize in wireless networks lll7h . we only consider bidirectional links in the 



secondary network when we define connectivity. As a consequence, when we determine whether 
there exists a communication link between two secondary users, we need to check the existence 
of spectrum opportunities in both directions. 

To summarize, under the disk signal propagation and interference model, there is a (bidirec- 
tional) link between A and B if and only if (CI) the distance between A and B is at most 
Tp, (C2) there exists a bidirectional spectrum opportunity between A and B, i.e., there are 
no primary transmitters within distance Ri of either A ox B and no primary receivers within 
distance r/ of either A ox B. 

2) Connectivity: We interpret the connectivity of the secondary network in the percolation 
sense: the secondary network is connected if there exists an infinite connected component a.s. 

Based on the above conditions (CI, C2) for the existence of a communication link, we can 
obtain an undirected random graph ^(A^, Xpr) corresponding to the secondary network, which is 
determined by three Poisson point processes: the secondary users with density A5, the primary 
transmitters with density Apt, and the primary receivers with density Apr (correlated to the 
process of the primary transmitters]^. See Fig. |4]for an illustration of ^(A^, Ap^)- 

The question we aim to answer in this paper is the connectivity of the secondary network, 
i.e., the percolation in ^(As, Apr). 

III. Analytical Characterizations of the Connectivity Region 

Given the transmission power and the interference tolerance of both the primary and the 
secondary users (i.e., Rp, Rj, Vp, and rj are fixed), the connectivity of the secondary network 

^The two Poisson point processes of the primary transmitters and receivers are essentially a snap shot of the realizations of the 
primary transmitters and receivers. In different time slots, different sets of primary users become active transmitters/receivers. 
Thus, even if a secondary user is isolated at one time due to the absence of spectrum opportunities, it may experience an 
opportunity at a different time and be connected to other secondary users. 
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Fig. 4. A realization of the heterogeneous network. The random graph Q{\s,^pt) consists of all the secondary nodes and 
all the bidirectional links denoted by solid lines. The solid circles with radii Ri denote the interference regions of the primary 
transmitters within which secondary users can not successfully receive, and the dashed circles with radii ri denote the required 
protection regions for the primary receivers within which the secondary users should refrain from transmitting. 

is determined by the density Xs of the secondary users and the density Apr of the primary trans- 
mitters. We thus introduce the concept of connectivity region C of a secondary network, which 
is defined as the set of density pairs (A5, Apr) under which the secondary network ^(A^, Apr) 
is connected (see Fig. [!])■ 

C = {{\s, Xpt) '■ ^(As, Apr) is connected.}. 

A. Basic Properties of the Connectivity Region 

We establish in Theorem [H below three basic properties of the connectivity region. 
Theorem 1 : Basic Properties of the Connectivity Region. 
Tl.l The connectivity region C is contiguous, that is, for any two points (A^i, Apn), (A5'2, Apr2) ^ 

C, there exists a continuous path in C connecting the two points. 
T1.2 The lower boundary of the connectivity region C is the As-axis. Let X*PTi^s) denote the 
upper boundary of the connectivity region C, i.e., 

X*prp{Xs)= sup{Apr : Q{Xs,Xpt) is connected.}, 
then we have that Apr(As') is monotonically increasing with A5. 
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T1.3 There exists either zero or one infinite connected component in ^(A^, Xpx) a.s. 

Proof: The proofs of Tl.l and T1.2 are based on the coupling argument, a technique 
frequently used in continuum percolation 11131 Section 2.2]. The proof of T1.3 is based on the 
ergodicity of the random model driven by the three Poisson point processes of the primary 
transmitters, the primary receivers, and the secondary users (the concept of ergodicity of a 
random model is reviewed in Sec. IV-A5I) . The details of the proofs are given in Sec. IV-B[ ■ 
Tl.l and T1.2 specify the basic structure of the connectivity region, as illustrated in Fig. [T] 
T1.3 implies the occurrence of a phase transition phenomenon, that is, there exists either a unique 
infinite connected component a.s. or no infinite connected component a.s. This uniqueness of the 
infinite connected component establishes the strong connectivity of the secondary network: once 
it is connected, it is strongly connected. It excludes the undesirable possibility of having more 
than one (maybe infinite) infinite connected component in the secondary network. We point out 
that such a property is not always present in wireless networks. Two examples where more than 
one infinite connected component exists in a homogeneous ad hoc network can be found in dlSJ. 



B. Critical Densities 

In this subsection, we study the critical density of the secondary users and the critical 
density Ap^ of the primary transmitters. Recall that 

A5 = infjAs : 3ApT > s.t. Q{\s,Xpt) is connected}, 
\*prp = sup{ApT : 3A5 < 00 s.t. Q{\s, Apt) is connected}. 

We have the following theorem. 

Theorem 2: Critical Densities. 
Given Rp, Rj, Vp, and rj, we have 

T2.1 \*g = Ac(rp), where Xc{rp) is the critical density for a homogeneous ad hoc network with 
transmission range {i.e., in the absence of the primary network). 

T2.2 X*prp < 4p^ax{'flV^}-r-^ ' ^'^^^^ '•^^ constant Ac(l) is the critical density for a homogeneous 
ad hoc network with a unit transmission range. 

Proof: The basic idea of the proof of T2.1 is to approximate the secondary network 
Q{Xs,Xpt) by a discrete edge-percolation model on the grid. This discretization technique is 
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often used to convert a continuum percolation model to a discrete site/edge percolation model 
(see, for example, [Il3l . Chapter 3], |40). The details of the proof are given in Sec. IV-C1[ 

The proof of T2.2 is based on the argument that if there is an infinite connected component in 
the secondary network, then an infinite vacant component must exist in the two Poisson Boolean 
models driven by the primary transmitters and the primary receivers, respectively. The key point 
is to carefully choose the radii of the two Poisson Boolean models in order to obtain a valid 



upper bound on Ap^. The details of the proof can be found in Sec. IV-C2I 




Fig. 5. A realization of the Poisson heterogeneous network when the percolation occurs (black stars denote primary transmitters, 
green plus signs denote primary receivers, red dots denote secondary users, and blue segments denote the bidirectional links 
between secondary users). We have removed secondary users who do not see opportunities for clarity. The simulation parameters 
are given by Apt ~ 10km~^, = 50m, Ri = 80m, As = 650km~^, Vp — 50m, ri = 80m, and the critical density in this 
case is Ac (50) ~ 576km~'^. 

Fig. [5] shows one realization of the Poisson heterogeneous network when A5 is slightly larger 
than Ac(rp) and Ap^ is small. At least one left-to-right (L-R) crossing and at least one top-to- 
bottom (T-B) crossing can be found in the square network. It is thus expected that these L-R and 
T-B crossings in finite square regions can form an infinite connected component in the whole 
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network on M^. If we slightly increase \pt, then we observe from Fig. [6] that the reduction in 
spectrum opportunities eliminates considerable communication links in the secondary network, 
creating several disjoint small components. 




Fig. 6. A realization of the Poisson heterogeneous network when the percolation does not occur (black stars denote 
primary transmitters, green plus signs denote primary receivers, red dots denote secondary users, and blue segments denote 
the bidirectional links between secondary users). We have removed secondary users who do not see opportunities for clarity. 
The simulation parameters are given by Apt = 20km"^ Rp = 50m, Ri = 80m, As = 650km"^, = 50m, n = 80m, and 
the critical density in this case is Ac (50) ~ 576km~'^. 

Fig. |7] shows a simulation example of the connectivity region, where the upper bound on the 
critical density X*prp of the primary transmitters given in T2.2 appears to be achievable. 

C. A Necessary Condition for Connectivity 

In this subsection, we establish a necessary condition for connectivity which is given in terms 
of the average conditional degree of a secondary user. This condition agrees with our intuition: 
the secondary network cannot be connected if the degree of every secondary user is small. 
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Fig. 7. Simulated connectivity regions when = 150m, rj = 240m, Rp = 100m, and Ri — 120m. The blue dashed line is 
the upper bound ^^^^^^i^l'i-^_^2 on the critical density Apy of primary transmitters given in T2.2. The area of the simulated 
heterogeneous network is 2000mx 2000m. For a fixed density As of the secondary users, the upper boundary Apy(As) is equal 
to the minimum density of the primary transmitters such that over all the 1000 realizations, the percentage of the ones in which 
there exists at least one L-R crossing is below 50%. The intuitive reason for choosing the existence of an L-R crossing as the 
criterion for connectivity is illustrated in Fig. I5I6I 

Let c?,rx/tx) denote the event that there exists primary receivers/transmitters within dis- 
tance d of a secondary user A. Let 1(^4, d, rx/tx) denote the complement of 1(^4, d, rx/tx). Since a 
secondary user is isolated if it does not see a spectrum opportunity, we focus on secondary users 
who experience spectrum opportunities and define the conditional average degree n of such a 
secondary user A as 

/i = E[deg{A)\ I{A,ri,rx) nl{A, Rj,tx)], (1) 

where deg{A) denotes the degree of A, rj the interference range of the secondary users, and 
Rj the interference range of the primary users. Notice that the degree of A is the number of 
secondary users within the transmission range of A and experiencing opportunities. We arrive 
at the following necessary condition for connectivity. 

Theorem 3: A necessary condition for the connectivity of Q{\s,\pt) is /i > 1, where fi is 
the conditional average degree of a secondary user defined in ([U). 
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Proof: The basic idea is to construct a branching process, where the conditional average 
degree /i is the average number of offspring. This branching process provides an upper bound on 
the number of secondary users in a connected component. If fi < 1, then the branching process 
is finite a.s. It thus follows that there is no infinite connected component a.s. in Q{\s,\pt)- 
Details can be found in Sec. IV-DI ■ 
To apply the necessary condition given in Theorem [3l the conditional average degree /i 
of a secondary user A needs to be evaluated based on the network parameters. Let 5 be a 
secondary user randomly and uniformly distributed within the transmission range Vp of A. Let 
g{\pT,rp,ri, Rp, Rj) denote the probability of a bidirectional opportunity between A and B 
conditioned on the event that A sees an opportunity. Based on the statistical equivalence and 
independence of different points in a Poisson point process, the conditional average degree fi of 
a secondary user A is given by this conditional probability g(-) of a bidirectional opportunity 
between A and a randomly chosen neighbor multiplied by the average number of neighbors of 
A, i.e., 

H = (XsTcrl) ■ g{XpT, r^, r/, Rp, Rj). (2) 

The detailed derivation for ^ and the expression for g{-) are given in Appendix A. It is also 
shown in Appendix A that g(-) is a strictly decreasing function of XpT- Thus g^^{-), the inverse 
of g{-) with respect to Xpt, is well-defined. 

Combining ^ with Theorem [3l we obtain an outer bound on the connectivity region. Specif- 
ically, let fi{Xs, Xpt) denote the conditional average degree of a secondary user in G{Xs, Xpt)- 
Then those density pairs (A^, Apt) satisfying fi{Xs, Xpt) < 1 are outside the connectivity region. 

Corollary 1: Given Rp, Rj, r^, and rj, an outer bound on the connectivity region C is given 

by 

where g^^{-) is the inverse of the conditional probability g{-) with respect to Apr- 

D. A Sufficient Condition for Connectivity 

In this subsection, we establish a sufficient condition for connectivity, which provides an inner 
bound on the connectivity region and a criterion for checking whether a secondary network is 
connected. 
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Fig. 8. An illustration of the dependent site-percolation model C with side length d (solid dots denote sites, solid lines denote 
edges connecting every two sites, and dashed lines denote the squared partition). 



The sufficient condition for connectivity is established by using the discretization technique. 
The continuum percolation model is mapped onto a dependent site-percolation model C in the 
following way. As illustrated in Fig. [8l we partition into (dashed) squares with side length 
d and locate a site at the center of each square. Sites whose associated dashed squares share at 
least one common point are considered connected (as illustrated by solid lines in Fig. [8]). Thus 
each site is connected to eight neighborj^ (see the eight neighbors Oi,...,Os of site O in Fig. [8]). 
Let Bq be the associated dashed square of O, then O is occupied if there exists in Bq at least 
one secondary user who sees an opportunity. 

Since the largest distance between two points in two neighboring dashed squares is 2\/2d, 
it follows that if we set d = then for every pair of secondary users in two neighboring 

^For the commonly used square site-percolation model, each site has four neighbors. The site-percolation model constructed 
here can provide a better inner bound. 
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dashed squares, they are within the transmission range of each other. Based on the definitions 
of occupied site in C and communication link in the secondary network, we conclude that the 
existence of an infinite occupied component (a connected component consisting of only occupied 
sites) in C implies the existence of an infinite connected component in the secondary network. 

Due to the fact that spectrum opportunities are spatially dependent, the state of one site is 
correlated with the states of its adjacent sites. Thus, the above site-percolation model £ is a 
dependent model. Define the dependence range k as the minimum distance such that the state 
of any two sites at distance d> k are independent, where the distance between two sites is the 
minimum number of neighboring sites that must be traversed from one site to the other. Then 
the dependence range of £ is given by 

"8max{i?, + ^,r, + ^}" 



k 



Tp 



1. (3) 



Let denote the upper critical probability of £ which is defined as the minimum occupied 
probability such that if the occupied probability p > p*, an infinite occupied component 
containing the origin exists in C with a positive probability (wpp.). Since the dependence range 
A; of £ is finite, it follows from Theorem 2.3.1 [12] that < 1- Now we present the sufficient 
condition for connectivity in the following theorem. 

Theorem 4: Let Pc denote the upper critical probability of the dependent site-percolation model 
L specified above. Define 

/(r,i?„r,)=2 rtMl||ll£)dt, (4) 

io TTi?^ 

where Si{t, Rp, rj) is the common area of two circles with radii Rp and rj and centered t apart. 
Then the secondary network is connected if 



1 — exp 



exp {-XpTiT [Rj + rj - I {Rj, Rp, rj)] } > Pc 



Proof: The proof is based on the ergodicity of the heterogeneous network model and 
its relation with the constructed dependent site-percolation model C. Details can be found in 
Sec. ICT ■ 
By applying a general upper boundon the upper critical probability for a site-percolation 
model with finite dependence range [12, Theorem 2.3.1], we arrive at the following corollary. 
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Corollary 2: A sufficient condition for the connectivity of Q{Xs, ^pt) is 



^PT < , 9 77-^ 7^ In 



1 — exp ' 



rc[Rj + r]~I{Rj,Rp,rj)] 1 _ ' 
where I{Ri, Rp, rj) is defined in dH) and k is the dependence range of the site-percolation model 
defined in (|3]). 

IV. Impact of Transmission Power: Proximity vs. Opportunity 

In this section, we study the impact of the secondary users' transmission power on the 
connectivity and the conditional average degree of the secondary network. As has been illustrated 
in Fig. m there exists a tradeoff between proximity and opportunity in designing the secondary 
users' transmission power for connectivity. Specifically, increasing the transmission power of the 
secondary users leads to a smaller critical density of the secondary users, but at the same 



time, a lower tolerance to the primary traffic load manifested by a smaller critical density Ap^ 
of the primary transmitters. 

A. Impact on the Conditional Average Degree 

As discussed in Sec. IIII-Cl the expression for the conditional average degree fi can be 
decomposed into the product of two terms: AsTrr^ and g{XpT, rp, r/, Rp, Rj). The first term is the 
average number of neighbors of a secondary user, which increases with the transmission power 
Ptx of the secondary users (i.e., enhanced proximity). The other term g{\pT, rp, rj, Rp, Rj) is the 
conditional probability of a bidirectional opportunity, which decreases with pt^ due to reduced 
spectrum opportunities. This tension between proximity and opportunity is illustrated in Fig. |9l 
where we observe that the impact of ptx on proximity dominates when ptx is small (yU increases 
with Ptx) while its impact on the occurrence of opportunities dominates when pt^ is large (/x 
decreases with pt^)- 

Corollary 3: Let pt^ be the transmission power of secondary users and p the conditional 
average degree defined in ([T]), then under the disk signal propagation and interference model we 
have 

^ = ((ptx)"^^°) as Ptx 00, 

''Here we use the Big O notation: f{x) — 0{g{x)) as a; ^ oo if and only if 3 M > 0, xq > such that < Af |ij(x)| 

for all X > xo. 
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where a is the path-loss exponent. 

Proof: We show this corollary by deriving an upper bound on the conditional average degree 
jjL. Details can be found in Appendix B. ■ 
For a homogeneous network, the average degree of a user is Xtit"^, which increases with ptx 
at rate {ptxf'^"' ■ In sharp contrast, this corollary tells us that for a heterogeneous network, when 
Ptx is large enough, the conditional average degree of a secondary user actually decreases 
with Ptx at least as fast as {ptxY'^^"' ■ 



8 




transmission range r^ 

Fig. 9. Conditional average degree of secondary users vs transmission range of secondary users (r^ cx {ptx)~ , where 
Ptx is the transmission power of secondary users and a is the path-loss exponent, and simulation parameters are given by 
Apt = 2.5km-^ Rp = 200m, Ri = 250m, As = 25km'^ n = rp/0.8). 

B. Impact on the Connectivity Region 

From the scaling relation of the critical density [13, Proposition 2.11], we know that in a 
homogeneous two-dimensional network, 

— 2 — — 

Ac(rp) = Ac(l) (rp) oc {ptx) " , 

where the constant Ac(l) is the critical density for a homogeneous ad hoc network with a 
unit transmission range. Thus, if each secondary user adopts a high transmission power, then 
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\c{rp) reduces. It follows from T2.1 that the critical density of secondary users to achieve 
connectivity reduces due to the enhanced proximity. 



On the other hand, from the upper bound on the critical density Ap^ of the primary transmitters 
given in T2.2, we have that 

A^ = O {{ptxy^^") as ptcc oo, 

where we have assumed that rp = Prj for some P G (0, 1) under the disk signal propagation 
and interference modej^. Thus, when the transmission power ptx of the secondary network is 



large enough, the critical density X*pj, of the primary transmitters decreases with ptx at least as 
fast as (ptx)"^''" due to reduced spectrum opportunities. 



C. Optimal Design of Transmission Power 

Due to the tension between proximity and opportunity, there does not exist a transmission 
power of the secondary users that leads to the "largest" connectivity region (largest in the sense 
that its connectivity region contains all regions achievable with any finite transmission power 
Ptx of the secondary users). Thus, the optimal design of pt^ depends on the operating point of 
the heterogeneous network. For instance, when a sparse secondary network is overlaid with a 
primary network with low traffic load, a large ptx may be desirable to achieve connectivity. The 
opposite holds when a dense secondary network is overlaid with a primary network with high 
traffic load. 

Focusing on a sufficiently dense secondary network, we address the design of its transmission 
power for the maximum tolerance to the primary traffic. Due to its tractability and achievability 
indicated by simulation examples (see Fig. |7]), the upper bound on the critical density \*prp of 
the primary transmitters given in T2.2 is used as the performance measure. 

Theorem 5: Let r/ and Rj denote the interference range of the secondary and the primary 
users, respectively. For a fixed Rj, the upper bound on \*prp given in T2.2 is maximized when 
the primary and secondary networks have matching interference ranges: = Rj. 



'since the minimum received signal power required for successful reception is, in general, higher than the maximum allowable 
received interference power , the transmission range rp is smaller than the interference range rj, i.e., P < 1. 
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Proof: Since under the disk signal propagation and interference model, r„ = (3ri for some 



/3 G (0, 1), the upper bound on \*pj, can be written as 

Ac(l) 



Ai^qhj for ^ 
(4^^ for r, > Rj. 



Then the above theorem can be readily shown by finding the maximal point for the two cases: 
Ti < Ri and rj > Rj. ■ 



An example of the upper bound on Ap^ is plotted as a function of r/ in Fig. \T0\ Notice that 
there is a distinct difference in the slope on the two sides of the optimal point. As a consequence, 
the operating region of rj < Rj is preferred over that of rj > Rj when the optimal point rj = Rj 
cannot be achieved. We point out that the desired operating region of rj < Rj is the typical 
case of a secondary network coexisting with a privileged primary network. 




Fig. 10. An example of the upper bound on Xpj, as a function of rj (Parameters are given by Ri — 120m, Vp = 0.625r/). 



V. Proofs 

In this section, we present proofs of the main results presented in Sec. UnHTVl We start with 
a brief overview of several basic results in percolation and ergodic theory that will be used in 
the proofs. 
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A. Percolation and Ergodic Theory 

1) Poisson Boolean Model: Poisson Boolean model is a common model in continuum perco- 
lation IibI. Often referred to as B{X, p, A), the model is specified by two elements: a Poisson 
point process X on with density A and a radius random variable p with a given distribution. 
Under this model, each point in X is the center of a circle in M'^ with a random radius distributed 
according to the distribution of p. Radii associated with different points are independent, and 
they are also independent of points in X. Under a Poisson Boolean model, the whole space is 
partitioned into two regions: the occupied region, which is the region covered by at least one 
ball, and the vacant region, which is the complement of the occupied region. We define occupied 
(vacant) components as those connected components in the occupied (vacant) region. 

Assume that nodes in a homogeneous ad hoc network form a Poisson point process with 
density A and their transmission range is r. It is easy to see that the connectivity of this network 
can be studied through examining the occupied connected components in the corresponding 
Poisson Boolean model B{X, r/2, A). 

2) Sharp Transition in Two Dimensions: Phase transition is a well-known phenomenon in 
percolation. For the Poisson Boolean model in two dimensions, this phenomenon appears more 
remarkable in the sense that the critical density for the a.s. existence of infinite occupied 
components is equal to that for the a.s. existence of infinite vacant components. Let Ac(2p) 
denote the critical density for the Poisson Boolean model B{X, p, A), then we have that 

□ when A < Ac(2p), there is no infinite occupied component a.s. and there is a unique infinite 
vacant component a.s.; 

□ when A > Ac(2p), there is a unique infinite occupied component a.s. and there is no infinite 
vacant component a.s. 



The exact value of Ac is not known. For a deterministic radius p, simulation results [|l9ll indicate 



that A, 



in m, 



'2p ) ^ 0.36p ^, while rigorous bounds 0.192p ^ < Ac(2p) < 0.843p ^ are provided 



3) Crossing Probabilities: A continuous curve in the occupied region is called an occupied 
path. An occupied path 7 is an occupied L-R crossing of the rectangle {0 < a; < /i} x {0 < y < 
I2} if 7 intersects with both the left and the right boundaries of the rectangle, i.e., 7 fl ({x = 

0} X {0 <y < I2}) 7^ 0, 7 n {{x = k} X {0 < y < I2}) ^ 0, and the segment between the two 
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intersecting points is fully contained in the rectangle (see Fig. [TTVa)). Similarly, we define an 
occupied T-B crossing by requiring that 7 intersects with the top and bottom boundaries of the 
rectangle (see Fig. [TlTb)). Let 



a((/i,/2), A, L-R) 
aiiluh), A, T-B) 



Pr{3 an occupied L-R crossing of [0, li] x [0, 12]}, 
Pr{3 an occupied T-B crossing of [0, /i] x [0,/2]}, 



denote the two crossing probabilities in the rectangle [0, h] x [0, 12]. Then for a Poisson Boolean 
model B{X, p, A) in two dimensions with a.s. bounded p, we have Corollary 4.1] that for 

any k > 1, 

' 1, ifA>A,(2p); 
0, if A < Ac(2p). 

Due to the symmetry of the Poisson Boolean model, similar results hold for the T-B crossing 
probability (T((n, A, T-B). 



lim (T((A;n, n), A, L-R) 



(5) 




Fig. 11. An illustration of the L-R crossing (a) and the T-B crossing (b) in a rectangle {0 < s < l\\ x {0 < y < l2\. 



4) Dependent Edge-Percolation Model: Let £ be a square lattice on with side length d 
(see Fig. [72)) . In an edge-percolation model, every site in L is occupied but every edge in C 
exists with some probability p. An existing edge is often referred to as an open edge, and an 
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edge that is not open is called closed. When the states (open/closed) of edges are correlated, we 
have a dependent edge percolation model. 



i 


€> 


®- 


e 


> 


© 

£+ 






' 1 

. \ 


I ; 
1 

\d 

----^ 

d 


\ < 

^ 


► 

d 

2 

e 


< 

d 
2 

\' 


t ( 

£ 

^ 


1 — 


4 


f % 

^ 


► < 

^ 


^ 

O 


> < 

^ 


» — 


1 


» < 


► < 




> < 


\ — 



Fig. 12. Part of the lattice L together with its dual L (solid dots and solid segments are sites and edges in L, and hollow 
dots and dashed segments are sites and edges in C^). The dual lattice £^ is the |)-shifted version of £, which is used in 
the proof of T2.1. Since distinct edges in L are crossed by distinct edges in and vice versa, there is a one-to-one mapping 
from the edges of L to the edges of . In this case, we claim an edge in being open if and only if its corresponding edge 
(i.e., the edge that it crosses) in L is open. 



Consider a special case of dependent edge-percolation model L where the state of an edge 
e is only correlated with its six adjacent edges (edges that share a common point with e). We 
have the following known result. 

Fact 1: y, Proposition 1] 
For any collection {cj}"^^ of n distinct edges in £, we have 

Pr{(Ci = 0) n (C2 = 0) n • ■ ■ n (a = o)} < gt, 

where Ci is the indicator of Cj being open,, and q = 1 — p h the probability of an edge being 
closed. 

5) Ergodic Theory: The study object of ergodic theory is the so-called measure-preserving 
(m.p.) dynamical system [Vt, 5", /i, T), which consists of a set Vt, a a-algebra 5" of measurable 
subsets of VL, a nonnegative measure fi on (Vt, and an invertible m.p. transformation T : Vt ^ 
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VL such that n{T~^F) = fi{F) for all F G 9^. A set F G 3^ is said to be T-invariant if T'^F = F. 
Obviously, all T-invariant sets in 5" form a a-algebra. 

An m.p. dynamical system (^7, 5", fi, T) is said to be ergodic if the a-algebra of T-invariant 
sets is trivial, i.e., for any invariant set, either it has measure zero or its complement has measure 
zero. Another property of the m.p. dynamical system that implies ergodicity is called mixing: 
an m.p. dynamical system (fi, 3", /i, T) is said to be mixing if for all E, F E 3^, fj.{T^E fl 
F) — n{E)n(F) ^ as n — > oo. For a m.p. dynamical system which is a product of two m.p. 
dynamical systems, we have the following classical result in ergodic theory. 



Fact 2: [|22 



Theorem 2.6.1] 

The product system of a mixing m.p. dynamical system and an ergodic m.p. dynamic system is 
ergodic, that is, for a mixing (fi, 3", yU, T) and an ergodic (\E', L, u, S), the product system 
(fi X \E', 3" X £, jj, X u, T X S) is ergodic, where 3" x £ is the cr-algebra on x ^ generated 
by subsets of the form F x L {F E 3^, L E and jj x u is the corresponding product measure. 

The concepts of ergodicity and mixing can also be defined for a random model under a 
probability space (l],?", yu), where the m.p. transformation T is replaced by a transformation 
group {Sx : X eM.'^ or Z'^} indexed by MJ^ or Z'^. For a point process model, the transformation 
Sx is usually to shift the realization lu E Vl hy x. A random model under a probability space 
(fi,?", /i) is said to be ergodic if there exists a transformation group {Sx : a; G M'^ or Z"'} 
that acts ergodically on (n,3^,fi). A transformation group {Sx ■ x G M*^ or Z'^} is said to 
act ergodically if the cr-algebra of events invariant under the whole group is trivial, i.e., any 
invariant event has measure either zero or one. Moreover, a random model under a probability 
space (fi, 5", /i) is said to be mixing if there exists a transformation group {Sx : x G M'^ or Z'^} 
such that for all E,F E J, we have ^.{SxE fl F) - j^{E)jj,{F) — > as \x\ oo. One direct 
consequence of an ergodic random model is presented as below. 

Fact 3: For an ergodic random model (Vl, 5", /i), if an event E E 3^ invariant under the whole 
transformation group {Sx : x G or Z'^} occurs wpp., i.e., n{E) > 0, then it occurs a.s., i.e., 
li{E) = 1. 



B. Proof of Theorem 1 

1) Proof of Tl.l: To prove Tl.l, it suffices to show that for any two given points {\si, \pTi) 
and (A5'2, ^pt2) in C, we can find a path in C that connects these two points. In particular, the 
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path we constructed is given by a horizontal segment and a vertical segment as shown in Fig.[T3l 
where we assume, without loss of generality, that Xsi < -^52- 



(A 



52, ApT2) 



(A5i,A 



PTl 



(A52, ApTi) 



(a) ApTi < ApT2 Ac 



A pt'^ 



(A51, ApTl; 

• 



(As'25 ApTi) 
<D 



(A52, ApT2) 



(b) Apri > ApT2 



Ac 



Fig. 13. The continuous path connecting the two points (Agi, ApTi) and (As2, ApT2) in the connectivity region C. 



Consider case (a) in Fig. [13] where Apti < Apt2- Case (b) can be proven similarly. First 
we show every point {Xs^^pti) (A51 < Ag < As2) on the horizontal segment belongs to C. 
Let A' = A5 — Asi. A Poisson point process X with density A5 is statistically equivalent to 
the superposition of a Poisson point process Xi with density A51 and an independent Poisson 
point process X' with density A'. It follows that any realization of the heterogeneous network 
with densities A5 and Ap^^i can be generated by adding more secondary nodes to a realization 
of the heterogeneous network with densities Asi and XpTi- Thus, the existence of an infinite 
connected component in ^(A^i, Apti) implies the existence of an infinite connected component 
in Q{\s-, Apri)- We thus have that (As, Apti) £ C for (A51 < A5 < A52)- 

Now we know that the two end points (A525APT1) and (A525APT2) of the vertical segment 
belong to C. For a point (As25ApT') (Apj-i < Apt < ApT2) on the vertical segment, let A' = 
ApT2 — Apt, then any realization of the heterogeneous network with densities As2 and Apt can 
be obtained by independently removing each primary transmitter-receiver pair with probability 
A'/ApT2 from a realization of the heterogeneous network with densities A52 and ApT2- It follows 
from the definition of communication link in the secondary network (see Sec. III-Bll) that the 
existence of an infinite connected component in ^(A525 ApT2) implies the existence of an infinite 
connected component in Q{\s-, ^pt)- Thus, we have (A52, Apt) G C (Apti < Apt < ApT2)- 

2) Proof of Theorem 1.2: Suppose that (A5, Apt) € C {XpT > 0), then by using the coupling 
argument for showing that the vertical segment belongs to C in the above proof of T 1.1, we 
conclude that (A5, 0) G C, i.e., the As-axis is the lower boundary of C. 
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Suppose that A52 > \si > 0. In order to prove the monotonicity of X*pj,{Xs) with A5 it suffices 
to show that V Apr > 0, if (A51, Xpr) G C then (A5'25 ^pt) ^ C. This is a direct consequence of 
the coupling argument for showing that the horizontal segment belongs to C in the above proof 
of Tl.l. 

3) Proof of Theorem 1.3: We first establish the ergodicity of the heterogeneous network 
model. 

Lemma 1: The heterogeneous network model is ergodic. 

Proof of Lemma [7]- The proof of this lemma is inspired by the proof of the ergodicity of 
Poisson Boolean model [13, Proposition 2.8]. The difficulty here is that for the heterogeneous 
network model, we have two correlated Poisson point processes: the primary transmitters and 
the primary receivers. The definition of the shift transformation for the primary network model 
is thus more complicated than the standard Poisson Boolean model with a deterministic radius 
p. To prove Lemma [H we first show the ergodicity of the primary network model, and then we 
show the mixing property of the secondary network model. Since the primary network model 
is independent of the secondary network model, it follows from Fact [2] that the heterogeneous 
network model is ergodic. 

Let B'^ denote the Borel a-algebra in W^, and N the set of all simple counting measureqj on 
S*^. Construct a cr-algebra !N for N generated by sets of the form 

{ne N : n{A) = k}, 

where A E 'B'^ and k is an integer. A point process X can now be defined as a measurable 
mapping from a probability space {Q, 5", P) into {N, !N) I21I . Chapter 7]. The measure fi on 
N induced by X is defined as fi{G) = P{X-\G)), for all G G N. 

In order to define the shift transformation on it is convenient to identify (fi, 5") with 
{N, Jsf). Let uj{A) denote the number of points in A G B"', V G and T^. be the shift 
according to a vector x G W^. Then induces a shift transformation Sx : i7 — ^7 through the 
equation for every A G B*^, 

iSxiu){A) = u;iT-'A). (6) 



*A simple counting measure on 2'' is an integer-valued measure for which the measures of bounded Borel sets are all finite 
and the measure of a point is at most 1. 
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Let {VLpT, 3^PT, Ppx) be the probability space of the Poisson point process Xpx for the 
primary transmitters with density Apy. Let Vlpp be the product space HneN Yizez^ ^Rv 
primary receivers, where Cr^ = {{x,y) : x'^ + y'^ < Rp}. Then we equip i^pR with the usual 
product cr-algebra and product measure Ppp with all marginal probability measure being fiu, 
where fiu is a uniform probability measure on Cr^^. Finally, we set Qp = ^Irt x i^pp and equip 
rip with the product measure Pp = Prt x Prr and the usual product cr-algebra. It follows 
that the primary network model is a measurable mapping from Vtp into Nrt xVLrr defined by 
{iOpT, curr) — * {Xpt{u}pt)-, ^^pr), where Nrt is specified in the definition of the point process. 

The positions of the primary transmitters corresponding to {ujrt, ujrr) E ^Irt x ^Irr are 
easily known from lurt- For the primary receivers, the positions are obtained as follows. Consider 
binary cubes 

2 

K{n, z) := JJ(^i2-", {zi + 1)2""] for all n G N and z e 1? . 

i=l 

For each primary transmitter xrt, there exists a unique smallest integer no = no(xpT) such that 
it is contained in a binary cube K{no, z(no, xrt)) which contains no other primary transmitters. 
The relative position of xrt^ receiver with respect to xrt is then given by ujRR{nQ, z{nQ, xrt))- 
Let ei, 62 denote the unit vectors in R^, then the translation Tg- : — > (i = 1, 2) defined 
by X ^ X + Cj induces a shift transformation f/g - on Vt rr through the equation 

{Ue,ujRR){n, z) = ujrrIu, z - 2"ei), for i = 1,2. 

Hence Tg. also induces a shift transformation Tg- on ^Ir = VIrt x ^^rr as follows: 

TeX^p) = {Se.UjRT, Ue.UjRR), fOT i = l,2, 



where Se^ is defined in Q. By using techniques similar to the proof of Boolean models |!l3, 
Proposition 2.8], we have that the m.p. dynamical system {^Ip, 3^p, Pp, TgJ is ergodic. 

Since the transmission range Vp of secondary users is fixed, the probability space of the 
secondary network model is the probability space (^5, 3^s, Ps) for the Poisson point process 
Xs of secondary users with density Xs- It follows from the proof of Poisson point processes jlA 
Proposition 2.6] that the m.p. dynamical system (Vis, 3^si Psi S^J is mixing. 

Since the primary network model is independent of the secondary network model, the sample 
space of the heterogeneous network model can be written as the product of Hp and ils, ^-c, 
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Vt = VLp^iVts- We equip Vt with product measure P = Pp x Pg and the usual product a-algebra. 
Similarly, the translation Tg. (i = 1, 2) induces a transformation T^. on = f2p x fi^, which is 
given by 

Then it follows from Fact [21 that the product m.p. dynamical system {Vt, 9^, P, TgJ is ergodic. 
Since the cr-algebra invariant under the transformation group {Tz '■ z E Z"^} is a subset of the a- 
algebra invariant under the transformation Tgi, we conclude that {T^ : z E Z,'^} acts ergodically, 
i.e., the heterogeneous network model is ergodic. ■ 
Based on Lemma [H we have the following lemma. 

Lemma 2: The number of infinite connected component in ^(A^, Xpr) is a constant a.s., and 
it can only take value from {0, 1, oo}. 

Proof of Lemma |2]- Let K denote the (random) number of infinite connected components 
in ^(A^, Apt), then since for all k >0, the event {K = k} is invariant under the group of shift 
transformations, it follows from Lemma [T] and Fact [3] that the event occurs with probability or 
1. Consequently, we have that K is an a.s. constant. Then it suffices to exclude the possibility 
of K > 2. This is shown by contradiction, that is, if there exist K > 2 infinite connected 
components, then they can be linked together as one connected component wpp. The proof is 
inspired by the proof of Proposition 3.3 in lll3n . and a major difference is that here we need to 
consider the impact of the primary network on the connectivity of the secondary network. 

Suppose that there are K > 2 infinite connected components a.s. If we remove all the 
secondary nodes centered inside a box B = [— n, n]^, then the resulting secondary network 
should contain at least K unbounded components a.s. Let, for A C M^, ^[A] denote the graph 
formed by secondary nodes in A. Given a box B and e > 0, consider the event 

E{B, e) := {d{U, B) < — e for any infinite connected component U in ^[-B*^]} . 

Partition the box B into squares with side length a > and let Sa = {Si, S^} denote 
the collection of all the squares which are adjacent to the boundary of B. Clearly, for a box B 
and e > 0, we can find a = a{B,e) G (0,rp/v^) and r] = T]{a) > such that for any point 
X ^ B with d{x,B) < Vp — e/2, there exists a square S = S{x) E Sa for which we have 
supj^gg d{x, y) < Tp — 7]. This means that, if we center in each square of Sa a secondary node 
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and there are neither primary transmitters nor primary receivers within a bigger box B = [—n — 
max{r/, Rj}, n+max{r/, -R/}]^, then every infinite component U in Q[B^] with d{U, B) < Vp — e 
is connected to some secondary node in Sa- 

Let E{a,'r]) be the event that each square in Sa contains at least one secondary node and 
E(B) the event that there are neither primary transmitters nor primary receivers within B. Since 
E(a, 7]) depends on the configuration of secondary nodes inside the box B, E{B, e) depends on 
the configuration of secondary nodes outside B and the configuration of primary nodes, based 
on the independence of the primary network and the secondary network, we have 

Fr{E{B , e) n E{a,7]) n E{B)) =Pr{E{B,e))Pr{E{a,r]))Pj:{E{B)\E{B,e)). 

If E{B,e), E(a,r]) and E(B) all occur wpp., then there is only one infinite connected 
componeni0 wpp. By using arguments similar to the proof for Proposition 3.3 in [13], we have 
that there exists a large enough box B and e > such that Fr(E(B, e)) > 0. Obviously, 
Pr{E{a,r])) > 0. Moreover, it is easy to see that P{E{B)\E{B,e)) > P{E{B)) > 0. ■ 

Now we have that the number K of infinite connected components is equal to zero, one or 
infinity a.s. To exclude the possibility of = oo, we can directly apply the proof of Poisson 
Boolean models [ Ijl, Theorem 3.6] here, which is based on several combinatorial results. The 
details are omitted. 



C. Proof of Theorem 2 

1) Proof of T2.1: To prove T2.1, it suffices to show that 

(a) for any A5 < \c{rp), the secondary network is not connected for any Apt > 0; 

(b) for any A5 > Ac(rp), there exists a \*pj.(\s) > such that V \pT < Apy(As'), the secondary 
network is connected. 

From Sec. IV-A[ we know that for a Poisson homogeneous ad hoc network with density A 
and transmission range r, the necessary and sufficient condition for connectivity is A > Ac(r). 
Since the existence of an infinite connected component in the secondary network implies the 
existence of an infinite connected component in the homogeneous ad hoc network with the same 
density and the same transmission range, by using a coupling argument, we conclude that when 

'since a < rp/\/E, every secondary node in a square of 5a is connected to those secondary nodes in the neighboring squares. 
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As < Ac(rp), there does not exist an infinite connected component a.s. in the secondary network 
for any Apr > 0. This proves part (a). 

The basic idea of the proof of part (b) is to approximate the secondary network Q{\s, ^pt) by a 
discrete dependent edge-percolation model on the grid. This discrete dependent edge-percolation 
model C is constructed in a way such that the existence of an infinite connected component in 
C implies the existence of an infinite connected component in G{\s, Apr)- 

Construct the square lattice C on M? with side length d (see Fig. [T2|) . Note that each site in 
C is virtual and is not related to any node either in the secondary network or in the primary 
network. Next we specify the conditions for an edge being open in C, which is the key to the 
mapping from G{\s, Apt) to C 

For each edge e in C, let (xg, Ve) denote the middle point of e. Then we introduce three 
random fields A^,, B^,, and Cg, all associated with the edge e in C, where Ce = v4ei?e is the 
indicator of the edge e being open, A^, represents the condition (CI) of the distance between two 
users for the existence of a communication link in the secondary network, and Be represents the 
condition (C2) of the spectrum opportunity. Specifically, consider the Poisson Boolean model 
B{Xs, rp/2, As) where Xs is the Poisson point process generated by secondary users, then for 
a horizontal edge e, Ae = 1 if the following two events (illustrated in Fig. [14)) are true: 

(i) there is an occupied L-R crossing of the rectangle [xe—3d/A, Xe+3d/A] x [ye—d/A, ye+d/A] 
in B{Xs, rp/2, As); 

(ii) there are two occupied T-B crossings of the square [xe~^d/A, Xe — d/A]x[ye — d/ A, y^+d/A] 
and the square [xe + d/A, Xe + 3d/A] x [y^, — d/A, ye + d/A] in B{Xs, r.p/2, As); 

and Ae = otherwise. For a vertical edge e, the definition of Ae is similar, where the horizontal 
and vertical coordinates are switched in the above two events. 

Next we define the random field Be. For an edge e in C, Be = I if Ae = 1 and the following 
two events are true: 

(i) there is no primary transmitter within distance Rj of any secondary node of the three 
crossings in the definition of Ae', 

(ii) there is no primary receiver within distance rj of any secondary node of the three crossings 
in the definition of Ae', 

and Be = otherwise. It follows from the definition of communication link in the secondary 
network (see Sec. III-B1|) that if Ae = 1 and Be = 1, then the three crossings in B{Xs, rp/2, As) 
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Fig. 14. A realization where Ae — 1 for the edge e (hollow points are sites in £ and solid segments are edges in £). 

are also three crossings in G{\s, Apt)- 

Let Ce = AeBe, then we claim that the edge e is open if Cg = 1, and e is closed if 
Ce = 0. We observe from Fig. [14] that whether the edge e is open is correlated with the states 
of the other edges. This model C thus is a dependent edge-percolation model. Furthermore, 
as shown in Fig. [151 if there exists an infinite open connected component in C, then those 
crossings associated with the edges in the infinite component in C comprise an infinite connected 
component in Q(\s, Xpr)- As a consequence, by considering the uniqueness of the infinite 
connected component in G{\s, ^pt), we only need to prove the following lemma in order to 
show T2.1. 

Lemma 3: Let C{0) denote the open connected component containing the origin O in £. 
Then given A5 > Ac(rp), 3 D > 0, Apj. > such that ford = D and any XpT < Ap^, we have 

Pr{|C(0)| = 00} > 0, 

where |C(0)| is the number of edges in C(0). 

Proof of Lemma \3} For an arbitrary edge e in £, let q = Pr{Ce = 0}, then we have 

q = Pr{(Ae = 0) U (Be = 0)} < Pr{A = 0} + Fr{B, = 0}. 
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Fig. 15. Percolation in C (thick segments are open edges in C and thin segments are closed edges in £, and blue curves are 
those crossings associated with the open edges). 



From the result on the crossing probabilities given in ([5]), we know that when A5 > Ac(rp), 

0. 

Thus when A5 > Kirp), Ve > 0, 3 D > such that Pr{/le = 0} < f . 

Given A^, = 1, let Srj be the area of the region covered by the circles with radii Rj centered 
at those secondary nodes in the three crossings, and Srj be the area of the region covered by 
the circles with radii rj centered at those secondary nodes in the three crossings. Then we have 

Pr{i?e = I y4e = 1} = Pr{3 some primary transmitter in Sr^} 

+Pr{3 some primary receiver in Srj}. 



PrjAe = 0} = Pr{at least one crossing does not exist} 



< 



a 



+ 




3d d 



d d 



1-a 



as d — > cx), i.e., limPr{Ae = 0} 
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Since < {f + 2Ri + rp) {I + 2Rj + rp) anas',, < (f + + r^) (f + 2rj + r^), it 
follows from the basic property of Poisson point processes that 



Pr{Be = \ Ae = l} < 1 -exp 



+ 1 — exp 



, 3(i 

-Apt ( y + 2i?7 + rp 
,'3d 

-XpT (— + 2rj + r 



d 



+ 2Ri + Tp 



- + 2r/ + 



Obviously, lim Pr{i?e = | = 1} = for fixed d. Thus if we choose d = D, then Ve > 0, 

A PT — ^0 

3 Apy > such that 



1} < - for all Apr < Ap^^- 



Pr{fie = Me 
It implies that when d = D, for all Ap^ < Ap^., 

Fr{B, = 0} = Ft{A, = 0} + Fr{B, = \ A, = l}Fr{A, = 1} 

1} 

2e 

< 



< Pr{Ae = 0} + Pr{5e = I 
2e 
3"' 

Thus for:d = D and all Apr < Apy, we have 



q < Pr{Ae = 0} + Pr{Se = 0} < e. (7) 

From Fig. [Ml we can see that if d > max{4i?7 + 2rp, Arj + 2rp}, then the state of edge e 
is only correlated with its six adjacent edges and it is independent of other edges. In this case, 
by using the 'Peierls argumenj^' [Q, Chapter 1], we can show that if the probability of an edge 
being closed q < then 

Pr{|C(0)| = oo} > 0. (8) 

The proof of the above statement follows the proof of Theorem 3 in [|31] except that the upper 
bound on the probability of n edges all being closed is replaced by the one given in Fact [T] 

Thus by combining ^ with (|7]), we conclude that given As > Ac(rp), 3 D > 0, Apy > 
such that for fixed d = max{D, ARj + 2rp, Arj + 2rp} and any Xpt < Ap^., 

Pr{|C(0)| = oo} > 0. 

Notice that Apj. depends on D which is chosen according to the crossing probability and is 
determined by As. As a consequence, Ap-p is a function of As, i-e., App = Apy(As). ■ 

*The essence of 'Peierls argument' is to make use of the one-to-one correspondence between a finite open component in 
lattice C containing the origin O and a closed circuit in the dual lattice of £ surrounding the origin O. 
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2) Proof of T2.2: From the conditions for the existence of a communication link in the 
secondary network specified in Sec. III-B 1 [ we know that for every secondary node in an infinite 
connected component, there can exist neither any primary transmitter within distance Ri of it 
nor any primary receiver within distance r/ of it. In other words, every secondary node in an 
infinite connected component must be located outside all the circles centered at the primary 
transmitters and the primary receivers with radii Ri and r/, respectively. Thus, if there is an 
infinite connected component in the secondary network, then an infinite vacant component must 
exist in the two Poisson Boolean models B{XpT, RpT, ^pt) and B{Xpji, Rpr, Xpt) driven 
by the primary transmitters and the primary receivers, respectively. Here Rpt and Rpji are some 
appropriate radii which will be specified later. A natural choice for Rpx is Ri, but if we consider 
the counterexample given in Fig. [161 then we can clearly see that even if there is an infinite path 
in the secondary network, no infinite vacant component exists in the Poisson Boolean model 
B{XpT, Ri, p\p) driven by the primary transmitters. Similarly, counterexamples can be easily 
constructed for choosing Rp^^ = rj. 




Fig. 16. A counterexample for choosing Rpt = Ri- AH the secondary nodes in the infinite path are located outside those 
circles centered at the primary transmitters with radii 7?/, which form a series of rings surrounding the origin O, and there is 
no infinite vacant component in the Poisson Boolean model B{XpT, Ri, \pt) driven by the primary transmitters. 

Suppose there is an infinite connected component in the secondary network. Then we can 
find a sequence of secondary users {Si, S2, S3, ■ ■ ■} such that they comprise an infinite path 
starting from 5*1 (see Fig. [TTI) . 

Assume that Si and Si+i (i > 1) are two adjacent secondary nodes in the above infinite path. 
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Fig. 17. An infinite patii in tiie secondary network. Tiie daslied segments form an inner bound on tiie infinite vacant component 
in tlie Poisson Boolean model driven by the primary receivers. 



Notice that the distance di^ j+i between Si and 5*4+1 satisfies di^ j+i <rp< rj, where the second 
inequality r-p < r/ follows from the fact that the minimum transmission power for successful 
reception is in general higher than the maximum allowable interference power. 

As we know, all the primary receivers must be outside the two circles with radii r/ cen- 
tered at Si and S'j+i, respectively, as shown in Fig. [TSl Given e > 0, consider the rectangle 
dj^^ ^ between Si and Si^i. By a simple computation in geometry, we have 

that the minimum distance from all the primary receivers to the rectangle is \j r'j — "^-j — e. As 
illustrated in Fig. [iTl it implies that there exists an infinite vacant component in the Poisson 
Boolean model B ^Xp/j, \Jrj — ^ — e, Xpr^ driven by the primary receiven^. By recalling 
the known results in Sec. IV-A21 we thus conclude that for all e > 0, 

APT< (2y/r|-r2/4-e)"'A,(l). 

Let e — > 0, then it yields 

^' I ' p 

The other term z Acfl) in the upper bound is obtained by applying the same argument 
to the Poisson Boolean model driven by the primary transmitters. 

D. Proof of Theorem 3 

Consider the connected component Ca containing an arbitrarily chosen secondary user A. 
Assuming that \Ca\ > 1, we construct a branching process as follows. Notice that if \Ca\ > 1 

'This technique used here can also be applied to the case when > rj, where only the minimum distance from all the 
primary receiver to the bar between 5"; and Si+i needs to be recomputed. 
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ri : 


Vi 

i+l '> C 















Fig. 18. One edge {Si, Si+i) in the infinite path. 



where \Ca\ is the number of users contained in Ca, then A must see the opportunity, i.e., 



I(yl,r/,rx) n I{A, Rj,tx) is true. Call A the initial point (or 0-th generation) of the branching 
process. Then the children of A {i.e., the 1st generation of the branching process) are secondary 
users which satisfy the following two conditions: 

(i) it is within distance of A, where Vp is the transmission range of secondary users; 

(ii) there exist neither any primary receiver within distance r/ of the secondary user nor any 
primary transmitter within distance Rj of the secondary user. 

The n-th {n > 2) generation of the branching process are obtained similarly, and they are 
connected to their parents in the {n — l)-th generation of the branching process via bidirectional 
links. Obviously, all the secondary users in Ca are counted in the constructed branching process 
model. But some of them may probably be counted more than once, since we do not exclude the 
previous n generations (including generation 0) when we consider the n-th generation. Thus, this 
branching process gives us an upper bound on the number of secondary users in Ca- It follows 
that if the branching process does not grow to infinity wpp., then there does not exist an infinite 
connected component a.s. in Q{\s, Xpr), due to the stationarity of the heterogeneous network 
model. Since the conditional average degree n is the average number of offspring for every 
generation, the necessary condition follows immediately from the classic theorem for branching 
processes [|l2l Theorem 2.1.1]. 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY, AUGUST, 2009. 



38 



E. Proof of Theorem 4 

From the construction of the dependent site-percolation model £, we know that the existence 
of an infinite occupied component in C implies the existence of an infinite connected component 
in Q{\s, ^pt)- Then in order to obtain a sufficient condition for the connectivity of the secondary 
network, it suffices to find a sufficient condition for the existence of an infinite occupied 
component in C. 

Let p be the probability that one site is occupied. Then based on the definition of the upper 
critical probability pc of C, we have that if p > pc, an infinite occupied component containing 
the origin exists in C wpp. It implies that if p > pc, there exists an infinite connected component 
in the secondary network wpp. Since the event that there exists an infinite connected component 
in the secondary network is invariant under the group of shift transformations, it follows from 
the ergodicity of the heterogeneous network model (see Lemma [B that if p > pc, there exists 
an infinite connected component in the secondary network a.s. 

Based on the definition of occupied site in C, we have 



In the last step, Pr |l {A, r/, rx) fl I {A, Rj, tx) | has been obtained by setting the distance d = 
in the expression for the probability of a unidirectional opportunity between two secondary users 
with distance d apart given in Proposition 1 in [14]. 



We have studied the connectivity of a large-scale ad hoc heterogeneous wireless network in 
terms of the occurrence of the percolation phenomenon. We have introduced the concept of 
connectivity region to specify the dependency of connectivity on the density of the secondary 
users and the traffic load of the primary users. We have shown several basic properties of 
the connectivity region: the contiguity of the region, the monotonicity of the boundary, and 
the uniqueness of the infinite connected component. We have analytically characterized the 
critical density of the secondary users and the critical density of the primary transmitters; they 
jointly specify the profile of the connectivity region. We have also established a necessary and a 




IiRi,Rp,rj)]}. 




VI. Conclusion and Future Directions 
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sufficient condition for connectivity, which give an outer and an inner bound, respectively, on the 
connectivity region. Furthermore, by examining the impacts of the secondary users' transmission 
power on the connectivity region and on the conditional average degree of a secondary user, 
we have demonstrated the tradeoff between proximity and spectrum opportunity. In establishing 
these results, we have used techniques and theories in continuum percolation, including the 
coupling argument, ergodic theory, the discretization technique, and the approximation using a 
branching process. 

To highlight unique design tradeoffs in heterogeneous networks, we have ignored the fading 
effect and the mutual interference between secondary users. If we take into account these factors, 
then the received signal to interference-plus-noise ratios at two secondary users will replace the 
distance between them in the condition (CI) for the existence of a communication link between 
them. This will result in a random connection model with correlated links, where the correlation 
between links is due to the mutual interference and the condition (C2) on the presence of the 
bidirectional opportunity. Although the connectivity region can still be defined in the same way, 
there will be another tradeoff between proximity and mutual interference besides the tradeoff 
between proximity and opportunity. The combination of these two tradeoffs will significantly 
complicate the characterization of the connectivity of the secondary network. We hope results 
obtained in this paper serve as a first step toward solving this more complex problem. 
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Appendix A: Expression for Conditional Average Degree 

The expression for the conditional average degree of a secondary user is presented in the 
following proposition. 

Proposition 1: Let A5 and Apt be the density of secondary users and primary transmitters, 
respectively. Let rj and Ri be the interference range of the secondary and primary users, respec- 
tively, and Tp and Rp the transmission range of the secondary and primary users, respectively. 
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Then the conditional average degree of a secondary user is given by 

/i = (XsTcrl) ■ g{XpT,rp,rj,Rp,Ri) 
rp 2t { r 

= AsvrrJ / - exp <^ - Apt i^ir] + R] + I{Rj, Rp, n)) - Si{t, rj, rj) - Si{t, Rj, Ri] 
Jo ''^p \ L 



where 



Si2{,r,9,Rp,t,ri] 
ttRI 

Su2{t,Rj ,Ri) 



rdrdO 



dt, (Al) 



T^Rl 

Si{t,ri,r2) the common area of two circles with radii ri and r2 and centered t apart (see 
Fig. [T9la)). and Su2{'t,ri,r2) is the union of two circles with radii ri and r2 and centered t 
apart (see Fig. [T9lb)). Sj2{r,9, Rp,t,rj) is the intersection area between one circle with radius 
Rp and the union of the two circles with both radii r/ (see Fig. [I2c)). For Si2{r,9,Rp,t,rj), 
the two identical circles are centered t apart, and the other circle is centered at (r, 6), where the 
middle point of the centers of the two identical circles is chosen to be the origin O. 

The expressions for /(i?/, -Rp, r/) and Si{t, ri, can be obtained in explicit form, which can 
be found in [14, Appendix A]. The expression for Si2{r, 9, Rp, t, r/) depends on the expression 



for the common area of three circles which is tedious and is given in [|24ll . By applying the basic 
property of the exponential function to (|AH) . we can easily show that g{-) is a strictly decreasing 
function of Apt. 

Proof: Let Kg (A) denote the event that there exist exactly k neighbors of a secondary user 
A. We thus have 



/i = E[deg{A)\I{A,rj,rx)nl{A,Rj,tx)] 



= ¥.KMdeg{A)\ I{A,ri,rx) nl{A, Ri,tx) nKs{A)]] 

oo f\ 2\^ 

= Y^^^xs^r^ y ^^^pJ E[deg{A)\I{A,ri,rx)nl{A,Ri,tx)nKsiA)]. 

k=0 

When A; = 0, it is obvious that deg{A) = 0. When A; > 0, let Bi be a neighbor of A, and Ipj 



an indicator function for Bi such that l^j = 1 if I{Bi,rj,rx) fl 1(5^, i?/, tx) occurs and Ipj = 
otherwise. Then by considering the statistical independence and equivalence of the k secondary 
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Fig. 19. An illustration of ri, r2) (the common area of two circles with radii ri and r2 and centered i apart), <S(72(t, fi, ''2) 
(the union area of two circles with radii ri and r2 and centered t apart), and Si2{r,9, Rp,t,ri) (the intersection area between 
one circle with radius Rp and the union of the two identical circles with radii rj). 



users, we have 



E[deg{A) \ I{A, n, rx) n Rj, tx) n Ks{A)] 

k 

^E[1b,| I{A,rj,rx)r]I{A,Ri,tx)] 



i=l 



= A;E[lsi I I(A, r/, rx) n I{A, Rj, tx)] 



kPr{I{Bi,ri,rx) nl{Bi,Ri,tx)\ I{A,ri,rx) nl{A, Ri,tx)} 



Pr{I(Bi,r7,rx) nl(5i,/?7,tx) nl(A,r7,rx) nl(A^/,tx)} 
Pr{I(^,r/,rx) n I(A^/,tx)} 



It follows that 



^Fr{I{Bi,ri,rx) r]I{Bi,Ri,tx) r]I{A,ri,rx) nl{A,Ri,tx)} 



Fr{I{A,ri,rx)nl{A,Ri,tx)} 



(A2) 



According to the definition of spectrum opportunity, Pr{I(A, r/, rx) n I(A,i?/,tx)} can be 
obtained by setting the distance o? = in the expression for the probability of a unidirectional 
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opportunity between two secondary users with distance d apart given in Proposition 1 in lll4ll : 



Pr{I(A, ri, rx) n Ri, tx)} = exp[-ApTvr(r; + R\ - I{Ri, Rp, r^))]. (A3) 



Next we derive the expression for the probability of a bidirectional opportunity, i.e., Pr{I(i?i, rj, rx)n 



Rj, tx) n I{A, rj, rx) fl 1(^4, Rj, tx)}, which depends on the location of Bi only through its 
distance to A. Since Bi is uniformly distributed within distance Vp of A, the density function of 
the distance t between Bi and A is given by for < t < Vp. In this case, the probability of 
a bidirectional opportunity can be written as 



Pr{I{Bi, rj, rx) nl{Bi,Ri,tx) nl{A,ri,rx) r]I{A,Ri,tx)} 

r'p 2t 

= / — Pr{I(5i,r7,rx) n I(5i,i?7,tx) nl(A,r-/,rx) nl(/l,i?7,tx)| d{Bi, A) = t}dt,{A4) 

Jo "^p 

where the integrand can be written as 



Pr{I(5i,r/,rx) nl(fii,i?7,tx) nl(A,rj,rx) nl(Ai?7,tx)| d{Bi,A) =t} 



Pr{I(5i,i?/,tx) n I(A^/,tx)| I(5i,r/,rx) nl(A,^/,rx) nrf(5i,A) =t] 



Pr{I(5i,r/,rx) nl(A,r/,rx)| d{Bi,A) = t]. (A5) 



Next, we compute the two probabilities in (IA5I) one by one. Since the primary receivers admit 
a Poisson point process with density XpT, we have 

Pr{I(5i,r,,rx)nl(A,r,,rx)| d{B^,A) = t} = exp[-ApT(27rrf - r,, r,))], (A6) 

where Sj{t,rj,rj) is the common area of two circles with both radii rj and centered t apart 
(see Fig. [Ha)). 

Let XpT denote the Poisson point process formed by primary transmitters. If we remove from 
XpT primary transmitters whose receivers are within distance r/ of Bi or A, then it follows 
from Coloring Theorem [15, Chapter 5] that all the remaining primary transmitters form another 
Poisson point process with density Apt 1 — '^'^^''^'f ifa"*''''^'' ' where S^ir, 0, Rp, t, rj) is the area 
of the circle with radius Rp and centered at (r, 9) intersecting the two circles with both radii rj 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY, AUGUST, 2009. 



43 



and centered t apart (see Fig. [I^c)). We thus have 



Pr{I(5i,i?7,tx) nl(A,/?7,tx)| I{Bi,ri,rx)r]I{A,rj,rx)nd{Bi,A) = t} 
exp < — Apr 



1 - 



Su2{t,Ri ,Ri) 



exp { -XpT 



2'kR] - Si{t,Ri,Ri) 



Si2{r,9,Rp,t,ri] 
nRl 



rdrd6' 



(A7) 



Su2it,Ri ,Ri) 

where Su2{t, Ri, Ri) is the union of two circles with both radii i?/ and centered t apart (see 

Fig.iiab)). 

Substitute dMl IMl) into (lA5l) . we have 



Pr{I(Si,r/,rx) nl(Si,i?/,tx) nl(A,r/,rx) nl(A^/,tx)| d{Bi,A) = t} 



exp s — A 



PT 



27,{rj + Rj) - Slit, rj, rj) - Si{t, Rj, Rj) 



Si2{r,9,Rp,t,ri) " 

— rdrdfc' 

TcRl 



(A8) 



Su2{t,Ri,Ri) 

The expression for the conditional average degree ji thus follows by plugging (|A8I) into (IA4I) 
and then dMl IMl) into dMl). ■ 



Appendix B: Proof of Corollary 3 

I — I 

From [14, Appendix A] and Fig. [T9lb. c), we know that when rj > Rp + Rj, 

I{Ri,Rp,ri) = Rl 
Si2{r,9,Rp,t,rj 



TiRl 



rdrd9 = Su2{t, Rj, Rj) = ^-nRj - Si{t, Rj, Rj). 



Su2ii,Ri ,Ri) 
Substitute dHU |B2l) into (|AT]) . we have 



2t 

IX = XsTirl I — exp[-ApT(7rrj - Si{t, ri, r/))]dt. 

n 

Plugging the expression for Sj{t,rj,rj) (14, Appendix A] into (IB3I ) yields 



^ = XsTlTp 



exp 



-\pT \ T^rj — 2rj arccos 



2ri 



+ t\lrj 



1 



dt. 



(Bl) 
(B2) 



(B3) 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY, AUGUST, 2009. 



44 



By applying the inequality arccos(x) <\ — x for < x < 1, we have 

2t 



'o "^p 

'•f 2t 



XsTT 



' p 

2 



< 



\2 2 \2 2 



■ exp(— ApT'i'"/)dt 

exp(-ApT/5r^) - -^exp(-ApT/??^f) 



A 



PT 



X2 



dt 



where we have assumed that Vp = Prj (0 < /? < 1) under the disk signal propagation and 
interference model. Since r/ oc (ptxY^", we arrive at Corollary [3l 
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